AdS/CFT duality at strong coupling 



(N 
O 

O 



> 
in 



- 1—1 

X 

5-H 



Matteo Beccaria and Carmine Ortix 

Dipartimento di Fisica, Universita' di Lecce, 
Via Arnesano, 73100 Lecce 
INFN, Sezione di Lecce 



O . 

• Abstract 



We study the strong coupling limit of AdS/CFT correspondence in 
the framework of a recently proposed fermionic formulation of the Bethe 
Ansatz equations governing the gauge theory anomalous dimensions. We 
provide examples of states that do not follow the Gubser-Klebanov-Polyakov 
0^ ' law at large 't Hooft coupling A, in contrast with recent results on the 

quantum string Bethe equations valid in that regime. This result indicates 
that the fermionic construction cannot be trusted at large A, although 
it remains an efficient tool to compute the weak coupling expansion of 
anomalous dimensions. 



q ! 1 Introduction 

\Q ' The AdS/CFT correspondence is a duality between string theory on AdS^ x S 5 

and Af = 4 super Yang-Mills [T] El El E El • If the number of colors N -> oo 
"p^ ' (the planar limit), the correspondence predicts a matching between a finite 

superconformal four dimensional theory and free string theory in a peculiar 
O !* geometry. Both sides of the correspondence are generically not solvable except 

in very special limits. In details, the duality map sends special massive string 
states to composite operators in the gauge theory. The mass of the string states 
equals the anomalous dimensions A in the gauge theory. The non triviality of 
the two sides of the correspondence implies the crucial fact that any test of 
AdS/CFT must explore the non perturbative behavior on at least one of the 
related theories. 

In this paper, we discuss the large 't Hooft coupling A limit where exact 
calculations can be performed in the string theory that can be treated in the 
supergravity semiclassical approximation. In particular, one should be able to 
recover the well known Gubser-Klebanov-Polyakov (GKP) result predicting the 
scaling law E ~ 2^/n\ 1 / i for the energy of level n massive string states as 
A — > oo 

This kind of investigation requires to determine the non perturbative behav- 
ior of anomalous dimensions in the gauge theory. This apparently out of reach 
task has a chance of being performed thanks to the conjectured integrability of 
J\f = 4 SYM [SUZIIH]- It is well known that there are Bethe Ansatz equations 
computing the correct perturbative expansion of the anomalous dimensions of 
operators with classical dimension L up to the L-th loop term (§]■ Technically, 
the dilatation operator is interpreted as the quantum Hamiltonian acting on 
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suitable subsets of gauge invariant operators closed under perturbative renor- 
malization. Each operator is mapped to the configuration of a spin chain. In 
the thermodynamical limit, this allows to derive closed exact expressions for 
the anomalous dimension of particular operators ^3^3- At fixed L, it seems 
hopeless to extend such technique up to the far strong coupling region. 

On the other hand, recently new Bethe Ansatz equations have been proposed 
based on a fermionic formulation of the dilatation operator in SYM The 
relevant lattice model is a twisted version of the Hubbard model (HM) . These 
equations have been conjectured to be free of the above wrapping problems. It 
is an open problem to determine whether they capture the strong coupling limit 
of the gauge theory. 

In this paper, we explore these equations and determine in particular the 
strong coupling expansion of two special states, the so-called antiferromagnetic 
state (AF) J0| and the folded-string one (FS) These are states placed in 
the opposite sides of the spectrum and are in a sense the simplest to analyze. 

Our results predict an asymptotic scaling A ~ yf~X for both states and ap- 
parently are in contrast with the GKP law. This fact opens an interesting 
discussion as we now explain. 

In principle, given a particular gauge invariant composite operator in SYM, 
it is non trivial to find the corresponding massive state in the string theory. If 
it is identified, at least in the semiclassical approximation, then it is possible to 
compute its energy at strong coupling and compare with the gauge theory result. 
In the case of the AF state, the asymptotic law A ~ y/X is conventional wisdom. 
In addition, and as a consequence, it has been proposed a particular semiclassical 
string configuration in the so-called slow string (SS) limit and conjectured to be 
dual to the AF state |HI . It has an energy scaling like y/X, but as we shall point 
out, with a different prefactor compared with the Hubbard model calculation. 
We believe that this is an indication that the proposed duality AF <-> SS is 
incorrect. Indeed, the recent work suggest that the y/X behavior of the 
AF state is totally fake and should not be trusted. This conclusion is based on 
a semi-analytic analysis of the quantum string Bethe Ansatz equations (SBA) 
valid by construction in the strong coupling limit 

In the case of the FS state, the scenario is somewhat different. There is no 
available solution of the SBA equations. On the other hand, the dual string 
state is identified quite reliably. As we show, the strong coupling behavior of 
the semiclassical energy reproduces in this case the GKP law. Thus, one is led 
to believe that the unavailable solution to the SBA equations should also exhibit 
a A 1 / 4 scaling, again in contrast with the fermion model calculation. 

In the following section we present the results obtained in the framework of 
the HM leaving to the concluding section a summary of their implications. 

2 The twisted Hubbard model 

As we discussed in the Introduction, the knowledge of the anomalous dimen- 
sions in the gauge theory is made possible by the correspondence between the 
dilatation operator in the J\f = 4 SYM and the Hamiltonian of a spin-chain 
model. In particular, in the SU(2) sector, the Hubbard model should predict 
at all loops and non-perturbatively anomalous dimensions of gauge operators of 
fixed length L A key point about the Hubbard model is that at strong 't 
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Hooft coupling A = gy M N, ^ reduces to a theory of free lattice fermions. How- 
ever one still needs to determine how the weakly coupled gauge theory states 
flow to lattice fermion states. In the following we will consider an half-filled 
Hubbard Hamiltonian : 

H = H + — — Hi 

47T 

where Hq contains the on-site Coulomb-type interaction: 

L 

H =L~Y^ ni,i n Ll , n al = c ai 
»=i 

while H\ is the free fermion part: 

/L-l \ 



* = E £4, 



i u <T,i+l + e C a,L C a,l 



<r=T4 \i=l 



The parameter (f> = n/2 arises from an Aharonov-Bohm flux which is needed to 
reproduce BDS results [§] • In addition, the introduction of the twisting phase 
makes the Hamiltonian invariant under the shift: 

Ca,j -> e^ /L Co-j+l, i = 1, • • • , L - 1, 

c CT)i - e**( i+1 )/ £ ^i, (1) 

and related transformation for annihilation operators. Information about com- 
posite operators anomalous dimensions can be derived from the analysis of the 
full Hubbard model spectrum. We will consider only states with total spin 
zero and invariant under the shift Eq.Q. This is a remarkable reduction of 
the Hilbert space. For example it has been shown that in a half-filled L = 8 
lattice the total dimension is reduced from 4900 to 226 17 . The coupling 
dependent spectrum flow has been evaluated numerically through direct diago- 
nalization method for lattices from L = 4 JHj to L = 8 ^3- Several crossings 
of the coupling dependent levels have been observed. This apparent violation 
of the Wigner-Von Neumann non-crossing rule follows directly from nontrivial 
coupling dependent conservation laws and is a characteristic signature of the 
one-dimensional Hubbard model quantum integrability |19| . 

Here a note of caution must be introduced. The fermion Hubbard model in- 
troduces states with double occupancy that do not have a direct correspondence 
with the gauge theory composite operators. The role of these extra states is still 
unclear, thus we will consider below the coupling dependent behavior of states 
that at A = reduces to states with no double occupancy, i.e. the states that 
have classical anomalous dimensions A(A = Q) = L (the perturbative multiplct 
in Ref 17.:). Nevertheless as the 't Hooft coupling increases, these states will 
mix with all the extra states with double occupancy. We will assume an opti- 
mistic attitude, waiting for a better understanding of the role of the Hubbard 
model. 

The direct diagonalization method is clearly turning off for large dimensions 
of the Hilbert space. Nevertheless it is possible to investigate longer composite 
operators using the integrability of the one-dimensional Hubbard model shown 
in Ref. ED]. 
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After working out the S-matrix one needs to diagonalize the multi-particle 
system by a nested Bethe Ansatz. The result of this procedure, generalized to 
the case with the twisting phase, yields the following Lieb-Wu equations for an 
half-filled lattice: 



L/2 



Lq n — 2ir I n + 2 tan 1 



2 ( u j ~ Tr sin (9« + 4>)) 



(2) 



L/2 L 

2tt Jk = 2 tan~ 1 (u^ — u j ) — 2 tan~ 1 

3=1 



m— 1 



2(wfc ^ 77^ sin(g TO + </>)) 



where 4> = n/(2L) and the /„, are the Bethe quantum numbers. The quan- 
tum numbers q n are the Bethe momenta of the particles while the Ui are the 
spin rapidities (Bethe roots) and describe the spin state. At A = the Lieb-Wu 
equations simplify since the equations for the momenta and the spin rapidities 
decouple and can be solved successively. Together with this decoupling, there 
is a decoupling of the wave function into a charge and a spin part that is deter- 
mined by the Bethe Ansatz of a Heisenberg spin chain. Once the spin rapidities 
of a particular state are determined in the Heisenberg model, it is possible to 
follow the evolution as the 't Hooft coupling increases, and finally to determine 
the strong-coupling behavior of any state. 



3 The AF operator 

At A = the state with highest anomalous dimensions corresponds to the ground 
state of the one-dimensional Heisenberg antiferromagnet. The Bethe-Ansatz 
solution in the thermodynamic limit has been found by Hulthen in 1938 |21j . 
Instead, we determine the Bethe roots for L finite and then follow the evolution 
as the 't Hooft coupling increases by solving step by step the Lieb-Wu equations. 
From the knowledge at every step of the Bethe momenta, we can determine the 
weak-strong coupling flow of the anomalous dimensions for the antiferromag- 
netic operator. In particular by determining the asymptotic free fermion state 
at strong coupling, we can provide a series expansion in inverse powers of A: 

Af F (A,L) =ao(L)V\ + ax(L) + a 2 (L)^= + ... (3) 

We find that the asymptotic free fermion state corresponds to the ground state 
of the Hubbard hopping term shown in Fig^ This is the state where all positive 
energy levels are doubly occupied and thus we have: 

\^ F ) = ii n °uio) 

n=l CT =U 

where a) apri is the canonical creation operator in the momentum space. Since 
the dispersion relation for the hopping term : 

/ 27m 7T \ 

£ „ = 2cos^— + -), 
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Figure 1: Asymptotic free fermion state for the AF operator at L = 20. 
we can easily obtain the function ao(L) of the expansion Eq.lPJ): 

(T\ 1 • -1 7F 

ao(L)= 27 Sm 2L 

In addition, we can compute analytically ai(L) from first-order perturbation 
theory. Expressing the Coulombian term Ho of the Hubbard hamiltonian in 
momentum space, we obtain: 



ai(L)= (tpQ F \L- - ^2 &{n- n +m- m')a\ n a ]:nl a\ m a^ ml \^ F ) 

It is easy to verify that the latter expression reduces to: 

ai (L) ee {^ F \L - I £ ^.r^A^a^n = L I (|) 2 = \h 

Finally the next-to next-to leading order coefficient can be put as: 

a 2 {L) — 2V2n Saf,l L 

with (!>af,l numerically evaluated for finite L £ 4N up to L = 32 (see FigEJ 
It is interesting to compare our expansion with the prediction of the AF state 
anomalous dimensions calculated via the BDS multi-loops Bethe Ansatz equa- 
tions. In this case a series expansion in inverse power of L is available. Thus: 

A^ S (X, L) = b (X)L + bi(X) + 6 2 (A)i + . . . (4) 

The leading order coefficient corresponds to the solution of the BDS Bcthc 
Ansatz in the thermodynamic limit L — > oo 10. It reads: 

^- 1 + ^{7x 
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Figure 2: Behavior of the next-to subleading term coefficient at strong coupling 
and finite L S 4N. 



where 



/(*) = 



dk J (fe) Ji(fc) 



fc 1 



and the Jj,(a;) are the Bessel functions of the first kind. In order to obtain the 
next-to and the next-to next-to leading order of the expansion Eq. J2J , one has 
to take into account finite size correction to the thermodynamic limit. This 
can be done by constructing a non-linear integral equation for the BDS Bcthc 
Ansatz |22|. The subleading order coefficient comes out: 



h=l 



dk 2-7T e 



ik xi, 



jh(g)_ 

VAfc cosh k/2 



where H and Xh represent respectively the number and the position of the Bcthc 
fake roots. In our case it turns out that H = since H = L — 2M and we are 
considering the case with M = L/2 magnons. Thus 

&i(A) = 

Contrary the next-to subleading order coefficient results: 



62(A) 



12 '0 



where I u {x) are the modified Bessel functions of the first kind. We can match the 
two expressions in the double limit A, L — > 00. Within the Hubbard formulation, 
the anomalous dimensions Eq. can be expanded at large L to obtain: 



A 



H 
AF 



L 



1 f 1 

24L + °U 



\/A+ -L + (Saf,, 



(5) 
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where the coefficient <5af,oo ~ 0.0240(1) has been found from a simple polyno- 
mial extrapolation of ^af.l- On the other hand, the BDS prediction can be 
computed in the strong coupling limit. The subleading order coefficient can be 
expanded asymptotically as |2H| : 

i N 

7T 4 * — ' 

m— 1 

where 

_ (2m - 1)(2 2 "'+ 1 - 1) [(2m - 3)!!] 3 ((2m + 1) ,, MI = 1 
23™-i(m-l)! ^ 2 ™+i ' 1 

and C is the Riemann zeta function. In addition, it is easy to verify that the 
asymptotic value of the next-to subleading order coefficient is given simply by: 

b 2 (A -> oo) = ^ 

Inserting the latter expansions in Eq.Q we find that the two Ansatz have a 
perfect matching in the terms cx L but a discrepancy of a factor 2 arises in the 
term \f\jL. This shows that only in the thermodynamic limit the two different 
approaches give the same results. The discrepancy in the finite size term may 
be due both to the uncontrolled effect of the wrapping interactions which start 
to contribute at the perturbation order o(A i ) in the BDS description and to 
the different order of limits A, L — ► oo. More interestingly the leading behavior 
at strong coupling in both descriptions of the anomalous dimensions scales as 
A 1 / 2 . It seems difficult to recover the generic GKP law 0: 

A~2^A 1/4 

It has been thus proposed to identify the dual string state of the antiferromag- 
netic state with slow-moving string solutions |14| that have the same scaling 
behavior. However such string states scale as: 

A - VAL 

This implies that the string solution has a numerical discrepancy l/7r 2 in the 
numerical prefactor with the dual gauge operator. This problem has been ear- 
lier attributed to the different order of limits since on the string side one first 
takes the large A limit and then consider L large while in the BDS gauge side 
description one first takes L large and then extrapolate the perturbative in A 
result at strong coupling. We have demonstrated that also in the Hubbard for- 
mulation this discrepancy survives although one consider the same limit order 
of the string theory and wrapping interactions do not contribute. Along a dif- 
ferent route, one can start from the classical string theory at large A and derive 
string Bethe Ansatz equations which are expected to match string calculations 
including leading quantum corrections. The latter approach has been recently 
studied in Ref . and the leading behavior at large A for the antiferromagnetic 
operator has been shown to reproduce the GKP law. 
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4 The FS state 



It is well known that Bethe solution of the Heisenberg antiferromagnet admits 
the existence of states with spin rapidities that acquire an imaginary part. In 
this case the Bethe wave function is made up of plane waves with complex 
wave vector. These states are usually named 'Bethe string solution' and can be 
identified with quantum Bloch wall state |24l 125) . In the AdS/CFT framework, 
the Bethe string state with all imaginary spin rapidities has the properties 
of the dual to semiclassical string considered in . We will instead consider 
the state with complex spin rapidities that corresponds to the state with lowest 
anomalous dimensions in the perturbative multiplet i.e. the lightest state. For 
this state, it has been shown that the Bethe roots (the spin rapidities) condense 
in the thermodynamic limit in two symmetrical curves in the complex plane 
1 1 1 1 and it is thus named the double contour solution. Its string state dual is 
nothing but the folded string state described in Ref.|13j 

As for the antiferromagnetic operator, we determine at g = the spin ra- 
pidities considering lattices with L = 12, 20 and determined step by step the 
evolution at strong coupling of Ui as well as of the Bethe momenta q n . For 
both lattices we find that in the limit g — > oo, the Bethe momenta flow to a 
state, shown in Fig.JJJ), where the positive energy levels with mode numbers 
n = k, L — k — 1 arc doubly occupied and the negative mirror levels are empties. 
The other positive levels are singly occupied as well as their negative mirror lev- 
els, thus they do not contribute to the leading term of the state energy. Since 
the results are completely identical for the two lattices we conjecture that for all 
L = 4(2fc + 1) the pattern will be identical. From the dispersion relation for the 
hopping term of the Hamiltonian the leading term of the anomalous dimensions 
at strong coupling results: 

A° FS (A,L) = ^cos^-^ (6) 

7T ZL 

The next-to leading term is computed from perturbation theory. It is easy to 
verify that the non-vanishing matrix elements of (tpQ s \Hq\?Pq S ) can be cast as: 

A Fs(9, L) = (*Po S \L -jYl a T,n a t>« a i,m a b m ~ \ A ^ ^ a \, m a iA^Q S ) 

n,m n,m 

(7) 

The first two terms have been computed for the AF operator. The third term, 
contrary is a spin-flip that gives non zero contribution only between singly oc- 
cupied mirror levels. Thus we have: 



Al s (g,L)=L-± 



(8) 



Finally the series expansion for the folded string anomalous dimensions can be 
put as: 

. , r , \/2 7T r- 3L 2 -2L + 8 . 2V2tt _ 

A FS (g,L) = — cos — VA + — + S xx £ + ••■ , 9 

7r ZL AL g 
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Figure 3: Free fermion state for the FS. The modes n = ±8 are empty. We 
show only a particular spin configuration for the singly occupied levels. 



Similarly to the antiferromagnetic state, the next-to next-to leading order co- 
efficient 5\ t L can be computed from the numerical exploration of the Lieb-Wu 
equation. We find: 



= 0.597(1) 
= 0.953(1) 



On the other hand, the energy of a folded string with angular momentum 
J = L/2 reads g?]: 



lim 



Afs(A'L 2 , L) 
L 



F(X') = -K(q) 



4qX 



-I 1/2 



E( q y 



where q = q(X') < 1 is the solution of 



4V 

^2" 



(K(q) ~ E{q)f E{qf 



(10) 



(11) 



and K(q), E(q) are standard complete elliptic integrals of the first and second 
kind. By expanding at large A', the second equations states: 

and reinserting in the expression for A we obtain: 



lim 

L — ►oc 



A FS (A'L 2 ,L) 



F(X') 



A'- 



V2 



A 'l/4 



Similarly to the AF state we find a natural strong coupling discrepancy between 
the gauge side and string side. The different scaling behavior in A suggests that 
it seems dangerous to rely on the gauge Bethe Ansatz equations to estimate the 
strong coupling limit of general states. 
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5 Conclusions 



To summarize, we have computed in the Hubbard model framework the strong 
coupling expansion of the anomalous dimensions of the AF and FS states. In 
both cases we find a leading term A ~ VA at large A. In view of the existing 
results discussed in the Introduction, we can draw several conclusions. 

First, in the case of the AF state, the Hubbard model is totally unable to 
predict reliably the strong coupling. The GKP law is violated whereas there are 
strong indications in support of it from the exact solution of the string Bethe 
Ansatz equations. 

Second, as a consequence, the duality AF <-> SS appears to be in quite bad 
shape. It is possible that for some gauge theory composite operator the scaling 
A ~ VA holds. However, this does not seem to be the case, at least for the AF 
state. 

Finally, in the case of the FS state, the conclusions are somewhat depending 
on the (not yet available) solution of the string Bethe Ansatz. If GKP scaling 
were found to hold, then we could conclude that even in this case the Hubbard 
model is not predicting the correct strong coupling. Of course, there is also the 
possibility of a violation of the GKP law. However, This option would rule out 
the duality with the folded string configuration. In our opinion, this possibility 
is disfavored due to the perturbative agreement (up to the 3 loop discrepancy) 
between the string calculation and the (weak) perturbative expansion of the 
gauge theory in the BMN limit 
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